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Abstract. We present a prospective work undertaken on Spectro-Polarimetric INterferometry (SPIN). Our theo- 
retical studies suggest that SPIN is a powerful tool for studying the mass loss from early type stars where strong 
Thomson scattering is present. Based on Monte Carlo simulations, we computed the expected SPIN signal for 
numerous hot star spectral types covering a broad range of geometries and optical depths. The SPIN technique is 
based on the detection and comparison of the fringe characteristics (complex visibility) between two perpendicular 
directions of polarization. The most obvious advantage is its ability to determine the polarization distribution 
in spherical winds for which no detection of polarization is achievable by classical techniques. In particular, we 
demonstrate that the SPIN technique is very sensitive to the /3 parameter from the so-called '/3 velocity law' for 
optically thin winds. Moreover, the location where the bulk of polarization is generated can be defined accurately. 
The required sensitivity for studying main sequence OB star winds is still very demanding (inferior to 0.5%), 
but the signal expected from denser winds or extended atmospheres is well within the capabilities of existing 
interferometers. The visibility curves obtained in two perpendicular polarizations for LBVs or WR stars can difi'er 
by more than 15%, and their corresponding limb-darkened radii obtained by the fit of these curves by more than 
35%. The signal expected from the extended circumstellar environment of Be stars and B[e] appears also to be 
easy to detect, relaxing the required instrumental accuracy to 1%. For these spectral types, the SPIN technique 
provide a good tool to extract the highly polarized and spatially confined envelope contribution from the bright 
star emission. 

It must be pointed out that the astrophysical environments investigated here offer a large panel of SPIN observing 
conditions in terms of geometry and polarization degree. The behavior of the SPIN observables can be transposed, 
at least qualitatively, to other astronomical objects for which important local polarization is foreseen. 
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1. Introduction 

Mass-loss is an intrinsic characteristic of hot stars 
which eject a strong wind during their whole short 
life. The light from the central star can be strongly 
polarized by its close circumstellar environment, es- 
sentially by Thomson scattering. The mass ejection is 
mainly driven by the pressure of the intense radiation 
field mediated by resonant scattering. The young tech- 
nique of optical interferometry has proven its efficiency 
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Luminous Blue Variable (LBV) P Cyg l|Vakili et al. 1997j) . 
Without spatial resolution, spectro-polarimetry represents 
one of the best suited techniques to study any departure 
from spherical symmetry of the mass-loss (see for instance 
Taylor et al. 1991 and Wood et al. 1997). The detection 
of a jet-like structure in the binary /? Lyrae with an in- 
terferometer (cf. Harmanec et al. 1996) and a spectropo- 
larimeter (cf. Hoffman et al. 1998) illustrates the comple- 
mentarity of both techniques. 

However, the interpretation remains limited by the av- 
eraging of the polarized information over the field of view 
since any observation of nearly symmetrical object pro- 
vides an almost undetectable signal. For instance, the pre- 
cision of current photo- and spectro-polarimetric observa- 
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tions is insufficient to test wind models with latitudinal 
dependance of the mass-loss rate in O supergiants, which 
predict a continuum polarization of only 0.1 per cent at 
most (Harries et al. 2002). 

Within this context, it appears very attractive to 
equip a long- baseline interferometer with a polari- 
metric mode in order to apply the so-called Spectro- 
Polarimetric INterferometry (SPIN) technique. Such at- 
tempts have been performed since the very beginning 
of interferometry. The unique Narrabri intensity in- 
terferometer was used with a polarimeter in 1974 to 
give an estimate of the polarization-dependent diam - 
eter change of (3 Orionis ijHanburv Brown et al. 197^ . 
but the signal-to-noise ratio (SNR) limitations were 
well above the expected signal. The experiment was 
repeat ed in 1981 with the I2T interferometer on a 
Lyrae ijVakih et al. 198l[) , and in 1997 with the GI2T o n 
the Be star 7 Cassiopeiae (jRousselet-Perraut et al. 19971) . 
These observations showed that instrumental po- 
larizati on has to be carefully studied and con- 
trolled I Rousselet-Perraut et al. 19971) . The first theoret- 



ical studies on SPIN were performed in the frame of the 
Narrabri Interferometer experiment by Sams & Johnston 
(1974). Cassinelh & Hoffman fl975') investigated the con- 
sequences of Thomson scattering around hot stars on the 
diameter measurements in linearly polarized light (with a 
single baseline). In the outer regions of the star, the light 
becomes polarized perpendicular in a direction parallel to 
the limb of the star. Integrated over the apparent disk this 
polarization cancels out. In contrast to this, an interfer- 
ometer in polarization mode can detect a signal due to its 
sensitivity to the polarized flux in a preferred direction. 
The star appears smaller in the plane of polarization par- 
allel to the baseline than in the plane perpendicular to it. 
An illustration of this effect can be seen in Fig. ^ that 
shows theoretical iso-intensity contours on the disk of the 
star for two orientations of polarization. Rousselet-Perraut 
(1998) performed a theoretical study of the SPIN observ- 
ables based on simple models with spherical and elliptical 
scattering environments. He also presents a methodology 
that is useful to interpret the results in the paper. 
We intend to present an updated overview of the SPIN 
capabilities using first, state-of-the-art Monte Carlo sim- 
ulations, and second a review of the on-going projects in 
the field of long-baseline optical interferometry (i.e. in- 
strumental facilities, signal performances to be obtained 
or expected in the near future). We develop a set of as- 
trophysical examples in the field of hot star winds which 
covers a broad range of wind geometries and intrinsic pa- 
rameters such as the wind density structure. In Sect. 13 the 
observables provided by long-baseline optical interferome- 
try are presented, together with the ones more specific for 
the signal study in polarized light provided by the SPIN 
technique. Sect. Ogives a brief description of the Monte 
Carlo code MC3D used and adapted for this purpose. In 
Sect. 01 we deal with spherical winds and perform numeri- 
cal tests for stars with different spectral types showing sig- 
nificant winds ranging from A supergiants to Wolf-Rayet 



stars. In Sect.[Sl we examine 2D geometries, ranging from 
anisotropic radiative winds to the disks of Be stars. We 
then discuss instrumental polarized devices foreseen and 
needed for such a technique in Sect. El Finally, we present 
the conclusions of this work. 

2. SPIN description 

In this section we describe the formalism applied in this 
article. At first, we briefly recall the interferometric ob- 
servables extracted from natural light. We restrict ourself 
to the case of a single interferometric baseline, i.e. with 
two telescopes. We adopt the formalism of Domiciano de 
Souza et al. (2002) and reproduce here the equations nec- 
essary for an introduction to natural light interferometry. 
In Sect. l2.2l this formalism is then extended to the case of 
polarized light. 



2.1. Natural light 

We consider a spherical star deflned by its hydrostatic ra- 
dius i?c, located at the center of the Cartesian coordinate 
system (x, y, z) shown in Fig. ^ The y axis is defined as 
the North-South celestial orientation and the x axis points 
towards the observer. 

Let us define the sky-projected monochromatic bright- 
ness distribution I\(y,z), hereafter called "natural light 
intensity map" . Interferometers measure the complex vis- 
ibility, which is proportional to the Fourier transform of 
I\{y,z). By denoting the Fourier transform of the inten- 
sity map hy I\{y, z) we can write the complex visibility in 
natural light as: 

VifyJ.A) - |V^(/„/.,A)|eW-/-^) = 

I\{v, 0) 

where fy and fz are the Fourier spatial frequencies associ- 
ated with the coordinates y and z. In long-baseline inter- 
ferometry the spatial frequencies are given by -Bproj , 
where XcS is the effective wavelength of the spectral band 
considered and Bproj is a vector representing the baseline 
of the interferometer projected onto the sky. The vector 
■Bproj defines the s direction, which forms an angle ^ with 
the y axis so that: 

-Bproj = (-Bproj cos C) ^ + (-Bproj sin ^) Z. (2) 

where y and z are unit vectors. 

We consider linear cuts along the Fourier plane corre- 
sponding to a given baseline direction J. We can define the 
new spatial frequency coordinates (u,v) for which -Bproj is 
parallel to the unit vector u. In that case the line integral 
(or strip intensity) of I\{s,p) over p for a given ^ can be 
written as: 



ds, 



The complex visibility is given by: 

F,(u,A) = |F,(u,A)|e*^^(-^) = ^. 

-^A,c(0) 



(3) 



(4) 
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Fig. 1. Adopted reference system. The figure represents 
isocontours on the apparent disk of a spherical wind as 
seen with a polarizer parallel (dashed line) and perpen- 
dicular (solid one) to the sky projected baseline Sproj- 
This axis forms an angle ^ with the sky coordinate system 
{x,y,z) and defines a new sky projected coordinate sys- 
tem {s,p) for which the s direction is parallel to Bproj- 
This (s, p) system is also the frame used in the polarization 
analysis. 

By varying the spatial frequency (baseline length and/or 
wavelength), we obtain the so-called visibility curve. Eqs. 
13 and 21 say that the interferometric information along 
-Bproj is identical to the one-dimensional Fourier trans- 
form of the curve resulting from the integration of the 
brightness distribution in the direction perpendicular (p) 
to this baseline. 



2.2. Polarized light 

SPIN allows to derive the geometry of the source as de- 
tected with the filtering view of the baseline and the po- 
larization. First, it must be stressed that contrary to clas- 
sical polarimetry, the polarizer direction is not fixed in the 
celestial North-South direction (y) but related to the base- 
line direction on the sky (s direction). This is due to the 
fact that the baseline is fixed to the ground, and not to a 
moveable mount as in the case of a monolithic telescopes. 
Thus it appears natural (and it is technically straightfor- 
ward) that the polarization analysis is using the baseline 
coordinate system (s, p) showed in Fig. ^ Throughout 
the entire article, two particular directions for this polar- 
izer are considered: polarizer in s direction (parallel to 
the baseline) and polarizer in p direction (perpendicular 
to baseline). 

Let /lin represent the polarized light contribution of the 
intensity map z). Our approach is adapted from the 
Stokes formalism to our time variable coordinate system of 
polarization analysis. We simplify the notation by writing 
/nat = Ix{y,z).We define 



-^s — -^nat -^lin 
-^p -^nat ~t~ -^lin ■ 



(6) 
(7) 



When the baseline direction is coincident with the celestial 
North-South direction, the polarization analysis system is 
coherent with the Stokes formalism: /lin = /q . 

The polarized visibility that is measured by a single 
baseline is defined by the Fourier transform of the strip 
intensity of the intensity maps modulated by the polar- 
ization (as defined in Sect. I2.1|l . When we observe with 
a polarimetric device, we record (simultaneously or not) 
three quantities: the visibility amplitudes curve |V| in nat- 
ural hght (iKiatl) and polarized light (|V^| and |V^|). These 
visibilities are related to the corresponding intensity maps 
Is and Ip by relations equivalent to Eqs. 01 and ^ These 
visibilities can then be related to a radius provided that 
a simple model of the object light distribution is defined 
(see Sect. I3.3|l . In order to study the polarized signal, we 
concentrate in this paper on a few observables chosen for 
their sensitivity and their simplicity for interpreting the 
geometry of the source. 

We define the polarized deviation curve AVp{f ) as the 
difference between the visibility curves in polarized light: 



A^p(/„/.,A) = \Vp{fyJzA)\ - mfyJz,X)\. 



(8) 



We can also define the degree of polarized visibility (fol- 
lowing the formalism from Rousselet-Perraut, 1997): 



\VpifyJ.,X)\~\VsifyJ^,X)\ ^ 



(9) 



The quantities AVp and Py are obtained at a given time 
for a given projected baseline. Three regions of the visibil- 
ity curves are of common interest in stellar interferome- 
try: the first lobe, the first minimum and the second lobe's 
maximum. 

The second lobe's maximum is very sensitive to the 
limb-darkening law of the star, and is consequently par- 
ticularly interesting for the study of diffuse light. However, 
precise observations in these high spatial frequencies re- 
quire long integration times in order to compensate the 
low fringe contrast. Moreover, as far as hot stars are con- 
cerned, the baselines required for studying the second lobe 
are generally larger than 200 m in the NIR band, except 
for the few closest stars. This is also true for the first mini- 
mum, commonly used for accurate radius determinations. 

Thus, it is more realistic to concentrate on the first 
lobe. We define the spatial frequency /max as the frequency 



where the SPIN signal AVp(/„ 



IS maximum. 



parallel -^perpendicular i 



and 



(5) 



The spatial frequency is expressed in units of the inverse 
stellar radius, 1/-R* and in the figures i?, = Re- As seen in 
the following sections, /max occurs at relatively low spa- 
tial frequencies, in the / — (0.2 — 0.6) l/Rc range, which 
relaxes the spatial resolution needed to perform SPIN ob- 
servations drastically. 
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3. The Monte Carlo code 

3.1. Presentation 

The simulation of visibilities and polarimetric observables 
are based on radiative transfer simulations performed with 
the Monte Carlo radiative transfer code MC3D (Wolf 
2003; see also Wolf et al. 1999, Wolf & Henning 2000). We 
assume a spherical, extended star which radiates isotropi- 
cally, i.e., the radiation characteristic at each point of the 
stellar surface follows the standard cosine law. The radi- 
ation field of the star is partitioned into "weighted pho- 
tons" each of which is characterized by its wavelength and 
Stokes parameters. The interaction of the stellar photons 
with the surrounding electron envelope is described by 
Thomson scattering. Due to (multiple) scattering events 
the polarization state of the initially unpolarized photons 
is modified. In order to derive spatially resolved images of 
the I, Q, and U Stokes vector components of the configura- 
tion, photons leaving the electron envelope are projected 
onto observing planes oriented perpendicular to the path 
of the photons. Since the optical depth in the electron 
envelope in some cases becomes <C 1, the enforced scat- 
tering concept introduced by Cashwell & Everett (1959) 
was applied in order to achieve a high signal-to-noisc ratio 
for the simulated images within a reasonable computing 
time. This concept has been used in particular in the study 
of the wind of ( Puppis presented in Sec. 14.11 and also in 
Sect, ma and Sect.lOI 

3.2. Limits 

In this study, an important limitation is that MC3D 
determines the polarization due to multiple photon 
scattering by electrons, but does not include the ef- 
fects of continuous hydrogen absorption and emission 
seen in disk-like circumstellar envelopes for example. 
Consequently, our modelling of a hydrogen disk in Sect.El 
results in an upper limit for the expected interferometric 
signal from an interferometer, especially for disk studies. 
This effect is discussed in Scct.|H| 



observed for each star. Diameters are usually derived by 
fitting the visibility of a well-defined artificial CLV like, 
e.g. a uniform disk (UD), a limb-darkened disc (LD) or a 
Gaussian intensity distribution to the observed visibility 
(cf. Jacob & Scholz 2002). The radius estimation based on 
the Uniform Disk (UD) assumption is not sufficient since 
we have chosen spectral types for which the diffuse light 
from Thomson scattering is important. A better way is to 
perform a fit of the visibilities based on the assumption 
that the emergent radiation follows a simple cosine limb- 
darkening law across the disk of the star (Limb-Darkened 
disk, LD): 



/ = 1 



(10) 



The two parameters of the fitted function are the LD an- 
gular diameter 8ld (©lds or 8ldp in polarized light) and 
the limb-darkening coefficient c (cg or Cp). These LD radii 
can be related by the canonical UD radius by (adapted 
from Sams & Johnston 1974): 



8 



UDs 



1 

1 - c/3 



(1/2) 



e 



LDs 



(11) 



In order to provide a clear view on the instrumental 
capabilities offered by contemporary interferometers, we 
have assumed a theoretical error curve based on a UD 
visibility curve estimation. Inside the first lobe a good an- 
alytical approximation for the visibility uncertainties aV 
due to the apparent star radius i?ap is given by (Vakili et 
al. 1997): 



aV 



\2h{z)\' 



(12) 



where J2 is the Bessel function of the first kind and second 
order. In the following, we assume that the interferometer 
is able to constrain i?ap with a 1% accuracy, and we build 
the related aV curve. This curve is then overplotted in 
each AVp diagram (in FigEl Figd FigEl Figlini FigCIl 
Fig ll3|l as a visual scale of the instrument accuracy for an 
illustration purpose. 



3.3. "Numerical" Diameters 

For each example treated in this paper, we perform a 
fit of the numerical visibility curves obtained with the 
code MC3D in order to derive the apparent diameter Oap, 
which is by definition larger than the diameter Qc de- 
fined with the hydrostatic radius Rc. The definition of 
the parameter Oap is not straightforward but depends on 
the model used for the fit. If the centre-to-limb variation 
(CLV) of intensity could be observed directly, an intensity 
radius could be defined in terms of the CLV shape and 
then related to a monochromatic optical-depth or filter 
radius via a model. In practice, reconstruction of the CLV 
from interferometric data is difficult with presently attain- 
able accuracies and the mean number of visibility points 



4. Spherical winds 

In this section, we consider spherical winds, for which the 
integrated polarized information averages out to a null 
value. In Tab.^ we present a list of typical early type stars 
for which the spatial characteristics (i.e. mainly the angu- 
lar diameter and the brightness) are well suited for inter- 
ferometric observations. Moreover, their winds are thick 
enough to expect a clear SPIN signal from their diffuse 
light. It must be pointed out that the angular diameters 
reported in this table are not the expected apparent an- 
gular diameters 8ap, but the core diameter or Oc used as 
input parameter for the radiative transfer simulation. The 
expected apparent angular diameter in natural light is in- 
creased by the wind as the star appears larger as a result 
of the diffused light from the halo. 
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Table 1. Some relevant parameters for the adopted models of spherical winds. The targets cover the range of spectral 
type for which strong local polarization are expected. Most of the parameters are adapted from Lamers & Cassinelli 
(1999), except for a Lyre and for WR 40 with parameters from Aufdenberg et al. 2002 and Herald et al. 2001. 









(" Puppis 


e Ori 


Deneb 


P Cyg 


WR 40 


Type 






04If 


BOIa 


A2Iae 


Blla/LBV 


WN8 


Distance 


D 


pc 


430 


410 


685 


1800 


2260 


Stellar radius 


Rc 


R© 


17 


35 


172 


76 


11 


Core angular diameter 


e 


mas 


0.35 


0.8 


2.35 


0.4 


0.04 


Stellar temperature 


T. 


K 


42000 


28000 


8875 


19300 


45000 


Mass loss rate 


M 


M© yr-i 


6 X 10"^ 


4 X 10"^ 


1 X 10"*^ 


1.5 X 10"'^ 


3 X 10"^ 


Terminal velocity 




km s~^ 


2200 


1500 


225 


210 


840 


Acceleration coefficient 


P 




1.0 


1.5 


3 


2.5 


1 


Optical depth 


To 




0.2 


0.17 


0.03 


1 


3.4 



For the wind velocity we apply the so-called /3 law: 



v{r) = Woo - ^) , 
where 

ro=i?c l-(^ 



(13) 



(14) 



The wind velocity is radial, accelerating to the value v^a- 
The stellar radius tq is defined as the hydrostatic radius 
i?c and in the following, wq is arbitrarily chosen to be 10 
km s^"'^. 

The mass loss rate is related to the density and the velocity 
via the following relation: 



M = Aiir^ p{r)v{r). 



(15) 



The local density is extracted from this relation for a given 
/3 law and mass-loss rate. 

C Puppis 

In order to illustrate the expected signal from a spherical 
wind, we modelled the star C Puppis (HD 66811), an early 
04If supergiant. Davis et al. (1970) have given an estima- 
tion of its apparent angular diameter based on a Uniform 
Disk fit: 8ap = 0.42 ± 0.03 mas. In the following, we per- 
form a study of the law parameters and then discuss the 
C Puppis 'standard' model. 

4.1.1. 'Standard' model 

The parameters of our 'standard model' are described in 
Tab. ^ Using this model, we find a maximum visibility 
deviation AV^nax=0.017 at the spatial frequency /max = 
0.47. Taking Agff —\[m\ and the angular diameter of C 
Puppis (Oap — 0.42 mas), this corresponds to a baseline 
of 300m, but only 180m at Aofr — 0.6/im. As we can see in 
Fig- El such a SPIN signal is close to the detection limit of 
an interferometer able to detect radius deviations of the 
order of 1%. 

Cassinelli & Hoffman (1975) also took C Puppis as ref- 
erence and gave the first estimate of the SPIN signal for 




Fig. 2. C Puppis fiux integrated in the baseline direction 
(strip intensity map) for polarizations s (solid line) and p 
(dotted line). Close to the photosphere, the wind density 
strongly decreases following a rapid acceleration modelled 
by the /3 law. 

a star with a wind. They used a two-component density 
model for which an extended atmosphere of total thickness 
Tc^lO is connected to a flow region with an optical thick- 
ness Tc=0.19, which is close to our standard parameters. 
However, their wind follows a different density law: 



(16) 



We have conducted a comparative study with the 
Cassinelli & Hoffman density law [n — —2) and the same 
To. Compared to our standard C, Puppis model (with the 
P law) , this new model presents a diffused light multiplied 
by a factor 2.5. For an equivalent optical depth, AV^ax 
is almost doubled, and /max = 0.18 only, to be compared 
with /max = 0.47 in case of our 'standard' model. The po- 
larization is generated much farther from the star with the 
power law model and the optimum baseline to detect the 
polarized signal taking Acg =1 /im (resp. 0.6/im) should 
be reduced to 120 m (resp. 70 m). 

By performing a least-square fit with a uniform disks 
to the simulated visibilities, Cassinelli & Hoffman expect 
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1.8 2.0 2.2 2.4 2.6 
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Fig. 5. Behavior of the maximum SPIN signal and the 
corresponding spatial frequency with parameter Voo- The 
dashed curve represents the course of the optical depth To. 



Spatial frequency (1 /R— *) 



Fig. 3. Visibility signal for both polarization, Vg in 
solid line, Vp in dotted line (top), and their difference 
A Vp (bottom) for the adopted parameter of C Puppis wind. 
The signal is small, but detectable if the interferometer 
sensitivity is such that an accuracy of 1% on radii mea- 
surements is possible (illustrated by the dashed curves, 
from Eq. [T^ . The discontinuity close to the spatial fre- 
quency / =1 appears because the first zero point is lo- 
cated at lower frequency for the polarization p intensity 
map, slightly more extended than the s one. 

an angular diameter ratio between two perpendicular po- 
larization directions of 7%, and Gap from a uniform disk 
fit of the star emission is 12% larger than the 0c. The 
results from the 'standard' model are smaller by a ratio 
of 2%, and 0ap=l.O70c. The contribution from the ex- 
tended atmosphere included in Cassinelli & Hoffman ap- 
pears to strongly favor the polarized signal compared to 
our Monte Carlo simulation of the wind component only. 
Their result has been confirmed later by Castor, Abbott 
& Klein (1975) who estimated that the scattering halo of 
Puppis increases Qc by 13% with a similar mass-loss 
rate {M — 6.6 x IO^^Mq/y). Nevertheless a refined study 
from Kudritzki et al.(1983) showed that the atmosphere of 
C Puppis cannot be considered as extended, which gives a 
lower limit Oc = 0.38 ± 0.03 mas, and a diameter increase 
of less than 10% in all cases and probably as low as a few 
percents, close to the result from the present 'standard' 
model. 

4.1.2. Mass-loss rate : M 

In this section we investigate how the mass-loss rate may 
affect the SPIN signal. As expected, we see in Fig.^that 
the SPIN signal AVijiax follows the mass-loss rate increase 
and the subsequent electron density increase almost lin- 
early. The wind from ^ Puppis is optically thin throughout 
the entire range of mass-loss rate encompassed. The spa- 
tial frequency /max at which this maximum AVmax can 
be detected is relatively stable. This means that the spa- 



tial location where the bulk of the polarization is gener- 
ated is relatively unaffected by a change of the mass-loss. 
However, the limb-darkening coefficients Cg/p, are very sen- 
sitive to the increase of diffused light since they are mul- 
tiplied by a factor larger than 3 as seen in Fig 01 This 
means that the relative light distribution, i.e. the balance 
of diffuse light near and far from the star evolves with the 
mass-loss. Slight multiple scattering effects are visible for 
high AI: AVmax evolves more slowly with M. 

4.1.3. Velocity : Voo 

Now M and (3 are kept fixed to their 'standard' values: 
M—6 X 10~^ and f3=l. In Eq. ^1 we see that increasing 
the terminal velocity Voo will decrease the local electron 
density and the scattering. The influence of a change in 
this parameter on the local electron density has a simi- 
lar impact as changing the mass-loss rate. The behavior 
of the SPIN signal follows the local electron density and 
the electron optical depth reported in Fig. [S] As men- 
tioned for the mass-loss rate, /max can be considered as 
constant. Nevertheless, changing Voo has a greater impact 
on AV;„ax than a change of the mass-loss rate: the slope 
of AVinax versus Voo in Fig.|31is twice as large as the slope 
versus M in Fig.^J The differences between the mass loss 
and velocity plots can be understood by the fact that the 
optical depth of the wind varies as Te = M /voo- Thus, Te 
and the visibility increase proportional to M and inversely 
with Vaa- Moreover the modification of Woo in the P law 
affects the local density close to the star (where the radia- 
tive field is strong) significantly by changing the value of 
ro . The scattering close to the star is increased by a larger 
amount by a change of Voo than by a change of M. 

4.1.4. Power-law index : P 

This parameter defines the density variation from the 
star to farther regions. Taking the parameters from the 
C Puppis standard model, we have varied /? between 0.8 
to 1.5. 
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Mass-loss rate (in 1 ^MQ.y_,) ^ Mass-loss rate (in 1 0"^fv1e.y_i) 

Fig. 4. Top: Behavior of the maximum SPIN signal AVinax and the corresponding spatial frequency /max with M . 
Middle and bottom: Result of fits of the polarized visibility curves with a 2 parameters limb-darkening law: LD angular 
diameters Glds (solid line) and 8ldp (dotted line) and LD coefficients Up (solid line) and Ug (lower line). In right the 
ratio 6ldp/0lds and Up/u^ are displayed. As the mass-loss increases, the contrast between Gldp and 6lds increases 
but Up/ Us evolves more slowly 



We see in Fig.lHl that /max and AVmax are strongly sen- 
sitive to p. This demonstrates that the parameter /max is 
very sensitive to the main location of the scattered light: 
The smaller /max, the flatter the density law, i.e. the bulk 
of the polarized emission is further away from the star. 
This behavior is particularly interesting since the /? pa- 
rameter is usually difficult to constrain by classical tech- 
niques using spectroscopic data such as spectrum fitting. 
A change of /3 has often a small influence on the spec- 
trum itself: the spatial extent of the line forming regions 
is affected but the lack of spatial information from spec- 
troscopy prevents from constraining efficiently this param- 
eter (see Hillier et al. 1998 for instance). Harries et al. 
(2002) pointed out recently the sensitivity of the polarized 
line profile morphology of O supergiants to the adopted 
velocity field. He suggests that it may be possible to use 
it as a diagnostic tool for the wind base kinematics. Such 
a diagnostic could also be performed by an interferom- 
eter with sufficient spectral resolution (see discussion in 
Sect©. 
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Fig. 6. Behavior of the maximum SPIN signal and the 
corresponding spatial frequency with parameter /?. The 
dashed curve represents the course of the optical depth 



4.2. € Oh 

The signal expected from e Ori is comparable to the one 
from C Puppis: Voo and M are similar, and thus Tc also. 
But this star is more interesting for interferometric obser- 
vations since the star is more extended than an 04If star 



with a similar expected SPIN signal. With Acff = 0.6/im 
(resp. Acff —1 /im, the baseline of half resolution (/ — 0.5) 
is only 70m (120m). This means that B supergiants offer 
a good compromise between polarized signal and spatial 
resolution. This is especially true for those which exhibit 
the strongest wind manifestations as LBVs (Sect.^3J. 
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Spotial frequency {1/R_*) 

Fig. 7. The expected signal polarization from the Deneb 
wind model is relatively weak, the AVp curve is below 
the 1% accuracy illustrated by the dashed curves defined 
from Eq.^] Nevertheless, a signal could be still detectable 
owing to the brightness and the angular diameter of this 
star. 



4.3. Deneb 

We performed simulations for Deneb (a Cyg ni) ba sed on 
the exhaustive work of Aufdenberg et al. l)2002|) . This 
A supergiant represents the "cool" detection hmit of 
Thomson scattering in the wind. The P law parameters 
are estimated from a fit of the nume rical electron den- 
sity law from Aufdenberg et al. l)2002l) using Eq. El The 
model parameters are displayed in Tab. ^ The optical 
depth due to electron scattering is about 7 times lower 
than for ^ Puppis, and the expected ratio between two 
perpendicular directions of polarization does not exceed 
0.5%. Nevertheless, we expect the results from our Monte 
Carlo simulation to be a lower limit for the expected SPIN 
signal produced in the extended atmosphere of this star 
(Sect. 13. 2() . and the radius ratio between perpendicular 
polarizations could be of the order of 1 %. The UD angu- 
lar diameter 9ap from Aufdenberg et al. l(2002^ of 2.4 mas 
is slightly larger than the core diameter expected from a 
non-extended atmosphere. However, it must be pointed 
out that Deneb's Ha profile exhibits a lack of the broad 
emission wing seen in the spectra of other supergiants, 
which are normally attributed to electron scattering. 
In natural light we expect that the apparent diameter of 
Deneb remains identical, whatever the baseline direction 
on the sky may be. The quasi-sphericity of the envelope is 
well established based on the absence of integrated polar- 
ization and detectable variation of radius with baseline for 
interferometric measurements (Aufdenberg et al. 2002). 

4.4. P Cyg 

LBVs are instable blue supergiants which cross a short- 
lived stage of instabilities with large mass loss rate in the 
HR diagram. The slow, dense wind from LBVs and in par- 
ticular P Cygni is at least optically thick due to Thomson 
scattering (tg = 1 for P Cygni), and multiple scattering 
occurs. The observed radius in natural light is very differ- 
ent from the hydrostatic one, defined by the basic P Cygni 
parameters. 



4.4.1. Smooth wind 

In Fig. |S1 the isocontours of the polarized intensity maps 
Is and /p are shown. They strongly depart from spherical 
symmetry, and this effect increases with the distance since 
the local polarization is larger in this case. The diffuse 
light is important and the angular diameter of the star in 
natural light represents 135% of the core angular diameter 
presented in Tab.^ Gap = 0.55 mas according to the pre- 
diction of the detailed P Cygni model of Najarro (2001). 
The SPIN signal is also large as shown in Fig El The op- 
timum baseline to detect the signal is 110 m (resp. 185 m) 
with Aeff = 0.6 /xm (resp. l^va). But the signal is still 
strong at shorter baselines: AVp=0.05 at a normalized 
spatial frequency /=0.15 (6% visibility difference), which 
corresponds to a 40 m baseline for Aog = 0.6/im and 70 m 
for Aoff — 1/im. The ratio between radii in perpendicular 
polarization directions reaches 1.24. The signal amplitude 
is well within the accuracy of optical interferometers which 
compensates the fact that the star is not well resolved. 

Such a strong SPIN signal is also expected for well- 
known LBVs such as AG Car, HR Car or HD316285. 
The last star has spatial parameters very close to those 
of P Cygni ones, i.e. an equivalent estimated distance, ra- 
dius or K magnitude, but differs by its huge mass-loss rate 
of the order of 2 x lO~^M0yr~^, more than 10 times larger 
than that of P Cygni (Hillier et al. 1998). Its wind is very 
thick (tc ~ 7) , and the ratio between radii in perpendicu- 
lar polarized lights reaches 1.37. Due to the high level of 
multiple scattering, the increase of mass-loss rate does not 
increase the SPIN signal compared to the case of P Cygni 
(AT/max=0.08). But the bulk of polarized light is generated 
far from the star and the weight of these extended polar- 
ized regions compared to the essentially unpolarized and 
heavily damped central star increases dramatically. The 
spatial frequency needed for resolving the polarized halo 
is consequently much lower compared to the one needed 
to resolve the central star, and the SPIN signal peaks at 
/max — 0.18. The apparent angular diameter in natural 
light based on the LD fit reaches 3.0x0c = 1.2 mas (2.5 for 
the UD angular diameter). It must be noted that the vis- 
ibility curve of HD 316285 is not adequately described by 
a Uniform or Limb-Darkened disk visibility fit. For these 
extreme spectral types fits based on the assumption of a 
gaussian distribution of light have to be preferred. 

4.4.2. Clumps 

The integrated polarization of P Cygni exhibits a strong 
variability (^^0.4%) on time-scales of days to weeks (Taylor 
et al. 1991), with no favored position angle which im- 
plies a quiescent state close to sphericity. This behaviour 
has been related to strong and localized eruptions and 
from an in-depth study of the spectopolarimetric vari- 
ability. Nordsieck et al. (2001) have strongly restricted 
the 'polarized' clump parameter space: position r < 2i?c, 
radius=0.1i?c, density contrast /9/po=20. These clumps 
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Fig. 8. Square root of polarized intensity maps contours 
for two perpendicular polarizers for the P Cygni smooth 
wind model p direction in dashed line and s direction in 
dashed-dotted line). The inner contour delimits the level 
0.81^^^, where /max is the maximum of the intensity /. It 
is close to the core angular diameter of 0.4 mas determined 
from P Cygni parameters defined in Tab ^ The outer 
contour delimits the level 0.3/,^„^. 
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Fig. 9. Visibihty Vg (solid Hne, up) and Vp (dotted Hue, 
up ) and the SPIN signal AVp (solid line, down) for the 
P Cygni smooth wind model. The SPIN signal is large and 
well above the 1% accuracy curves (dashed lines, down). 
Ay,nax=0.087 at /max = 0.40, which corresponds to a rel- 
ative signal Pv larger than 12%. 



could be related to those detected further out (about 0.5 
arcsec) with Adaptive Optics by Chesneau et al. 2000. 

Previous interferometric studies have demonstrated 
that the modulus of the visibility (as used in the paper) 
is not very sensitive to small scale asymmetries of the ob- 
ject (Vakili et al. 1997). However, the photocenter of the 
emission in natural and polarized light is no longer cen- 
tered on the star and consequently affects the phase of 
the fringes. The expected signal should be faint but a dif- 
ferential study between two polarizations or between line 
and continuum should increase the accuracy. Such a study 
for photometry and polarimetry, similar to the one carried 
out by Rodrigues & Magalhaes (2000), is not in the scope 
of this paper and deserves more extensive investigation. 

4.5. WR 40 

The Wolf-Rayet (WR) evolutionary stage is characterized 
by a strong mass-loss rate and a very dense and optically 
thick wind. These stars also exhibit a much smaller hy- 
drostatic radius, and a faster terminal velocity compared 
to the LBV stage. 

Most of the brightest WR stars are located in distances 
larger than 1 kpc and their core radii are generally smaller 
than 4 Rq which implies core angular diameters Qc lower 
than 0.04 mas. Nevertheless, the expected apparent diame- 
ter can be much larger because the wind is optically thick 
far from the star. Among WR subtypes, the WN8 star 
offers a good compromise between luminosity and core ra- 
dius extent, which reach 10-15 Rq, and they can therefore 
be more easily detected with long-baseline interferometers. 



The parameters of WR40 are extracted from the ded- 
icated study of Herald, Hillier & Schulte-Ladbeck (2001), 
and the SPIN signal presented in FigHni As for HD316285 
AVmax saturates to the P Cygni value, but the shape of 
AVp is characteristic of a strongly optically thick wind 
(see also Cassinelli & Hoffman 1975): /max=0.26, and the 
visibility curves can no longer be fitted by UD or LD disk 
visibility laws. 

The WR wind parameters are well constrained with 
current line-blanketed non-LTE models atmospheres. 
Nevertheless, the constrains on the /3 parameter are still 
weak. Herald et al. (2001) investigated (3 over a range of 
0.5-2 and they did not detected any significant impact on 
their spectra. By varying /3, we notice that the SPIN signal 
is also no longer sensitive to this parameter. The optically 
thick zone is so large that most of the wind acceleration 
is embedded in it and the polarization comes essentially 
from regions close to the terminal velocity. The study of 
the /3 law can only be conducted for optically thin winds. 



5. Disk geometries 

In this section, we investigate objects for which the close 
environment can no longer be considered as spherically 
symmetric and presents a 2D structure created by a co- 
latitude dependance of the mass-loss rate. The generated 
structure can be an extended compressed equatorial region 
like for B[e] stars or even a disk for Be stars for instance. 
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Fig. 11. Model of Tau for 3 inclinations corresponding to 3 disk inclinations (columns). The upper subpanels show 
the intensity maps on the sky with the overplotted polarization P. The thick bar represents a linear polarization of 
100%. The lower panels present the overplotted visibility curves for Is and Ip maps (s polarization in solid line and 
p polarization in dotted line), and their difference AVp for a vertical (z, middle panels) and horizontal (y, bottom 
panels) baseline orientation. For each of the AVp curves corresponds a couple of 1% sensitivity curves (dashed) defined 
in Eq. ^|to illustrate the strength of the signal. 



5.1. Be stars 

Be stars are hot and fast rotating stars surrounded by an 
extended circumstellar hydrogen envelope. They manifest 
the so-called 'Be Phenomenon' characterized by Balmer 
lines in emission and infrared excess. One of the challeng- 
ing questions on Be stars is the geometry of their disk, and 
in particular their opening angle, about which there is still 
an active debate. Most authors have considered geomet- 
rically thin disks (half opening angle of 2-5°). The very 
narrow disks considered by Wood et al. (1997) were those 
predicted by the Wind-Compressed Disks (WCD) theory 



of Bjorkman and Cassinelli (1993). Furthermore, interfer- 
ometric obse rvations have given uppe r limits of approxi- 
mately 20° llOuirrenbach et al. 1997t) . However, the hy- 
pothesis of such narrow disks faces several problems, and 
the current set of observations does not provide a unique 
interpretation on the circumstellar geometry (Yudin ct al. 
1998). 

The model from Waters l)l986|) has been successfully used 
to explain the near and far IR observations and is coherent 
also with polarization data (Cote & Waters 1996, Waters 
& Marlborough 1992). They model the disk as an equa- 
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Fig. 10. Visibility Vs (solid line, up) and Vp (dotted line, 
up ) and the SPIN signal AVp (solid line, down) for the 
WR 40 model. AVp is well above the curves of 1% ac- 
curacy (dashed lines, down). The second lobe of the vis- 
ibility function almost disappears by apodisation for the 
WR model because the central star is no longer detectable 
within the envelope. 



torial cone with a density law described in Eq. 1161 with 
a density gradient n ~ 2 — 3.5, the density of the disk 
po ^ 10^^ — 10^^ cm~^, the disk radius i?d, the viewing 
angle i, and the disk half-aperture 9 as main parameters. 
We performed simulations of the polarized emission of Be 
star disks. These simulations have been tested using the 
work from Wood et al. (1996, Fig. 112(1 . We obtain very sim- 
ilar results for the integrated polarization, (with less than 
10% deviation) taken into account that our density law 
is not exactly similar to that used by Wood et al. which 
allows a colatitude dependency of the density. However, 
their polar-to-equatorial density ratio of 1 : 10"^ is very 
strong, so that we can consider the two models being very 
close. In particular, we used as template the star Tau 
(Bl IVe-sh star), which is a well studied case of an edge- 
on Be star (i ~ 82°; Wood et al. 1997). We used a set 
of parameters close to those of Wood et al. (1997), i.e., 
i?* = 6i?0, T, = 2 X 10^ K for the central star and a dis- 
tance of 120pc. For the disk simulation a very thin disk of 
half-opening angle of 3° and density law which defines an 
optical depth in the disk plane of about — 3 (density 
gradient n = — 3) are used. We recall that in our study the 
only emission process considered is Thomson diffusion. 

In Fig. llll we present the expected interferometric sig- 
nal for C Tau for different inclinations. In the pole-on case, 
we clearly see that the visibility deviation curves AVp{f) 
are identical but inverted between both baseline (the base- 
line is aligned and perpendicular to the polarizer). This 
configuration is the most favorable for the local polariza- 



tion: the bulk of electrons is perpendicular to the line of 
sight and AT4iax reaches 0.03. For i — 45°, the local po- 
larization in the vertical direction is much weaker. This is 
due to the thinness of the disk: for most of the electrons, 
the polarization efficiency is only 30% at the diffusion an- 
gle X =45° (P = (1 - cos(x)^)/(l + cos(x)^)). The SPIN 
signal is almost undetectable but for the perpendicular 
baseline, the signal is much stronger {AVma^ = 0.045). We 
see in this example that the ratio of AVp between two per- 
pendicular baselines provides valuable information on the 
system inclination i on the sky and about the aperture of 
the disk. This complements the natural light information, 
i.e. the ratio of the radii in two perpendicular directions 
which provides information on the projected 2D intensity 
map but no indication on the real 3D structure of the ob- 
ject. 

In the equator-on configuration, an unexpected large sig- 
nal (|AI/max| = 0.02) is visible at a high spatial frequency 
(/max = 0.83) with a vertical baseline (i.e. perpendicular 
to the disk). This effect is related to the increasing disk 
vertical extension in the external regions and is strength- 
ened by the absorption of the unpolarized star light in the 
line of sight. With an horizontal baseline, /max — 0.37, 
i.e. lower than in the previous case, and AVmax — 0.03. 
This is also due to the optical thickness of the disk, which 
prevents the observation of high polarization close to the 
line-of-sight to the star. The bulk of polarization is thus 
located further out, and /max is decreased compared to 
other inclinations. 

Fig. 1121 illustrates the differences and the complemen- 
tarity between the polarimetric and interferometric ob- 
servables. We reproduce quantitat ively the integrated po- 
larization curve from Wood et al. l)l996|) . In edge-on view, 
there is a strong decrease of the local polarization due to 
the highly optically thick line of sight (multiple scatter- 
ing). Nevertheless the integrated polarization is compen- 
sated by the higher asymmetry of the system (right panel 
in the figure). In pole-on view the integrated polarization 
is null but the SPIN signal is still detectable and even 
slightly larger than in the equator-on case. For integrated 
polarization and SPIN signals the maximum occurs at the 
same angle i ~ 65°. The interferometer is very sensitive 
to the large projected emitting surface of the pole-on view 
which compensates the lower local polarization. 

5.2. B[e] stars 

B[e] stars are hot supergiants showing an important ex- 
cess in the infrared due to the presence of hot circum- 
stellar dust. These stars exhibit also the so-called 'Be 
Phenomenon', but also show forbidden lines in their spec- 
trum. Zickgraf et al. (1985) proposed a model for the LMC 
B[e] supergiant R126 consisting of a fast wind in the po- 
lar regions and a dense and slow wind in the equatorial 
region where the dust is formed. In contrast to Be stars, 
the circumstellar geometry of B[e] stars is rather more of 
an open question. Moreover, within the sample of clas- 
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Fig. 12. Left: Polarimetric and SPIN signal for various inclinations. Middle: Expected SPIN signal, represented by 
the maximum of the deviation curve (as shown in Fig|3J). Right: Ratio of uniform diameters (UD) fit between two 
perpendicular polarization direction. The upper (lower) curve describes the ratio for an horizontal (vertical) baseline. 
Due to the spherical symmetry, the pole-on signal is just inverted between the baseline. At higher inclinations rn 
increases, but the vertical polarized signal disappears and ry reaches almost 1. 



sified galactic B[e] stars which could be resolved by an 
interferometer, one can find young stellar objects such as 
extreme Herbig Be stars (HaebeB[e]), together with super- 
giant stars (sgB[e], see Lamers et al. 1998). The distance 
estimations and hence their luminosity and angular diam- 
eter are poorly constrained, which render the classification 
problem. As revealed from their emission lines and near 
IR (J band) excess, the envelopes of sgB[e] stars provide 
an ample opportunity for scattering of radiation from the 
central star by free electrons. Moreover, dust is evidenced 
by their IR excess (Melgarejo et al. 2001). It is possible but 
sometimes difficult to discriminate the electron scattering 
and the dust scattering regions with spectropolarimetry 
alone (Oudmaijer & Drew 1999) and our SPIN technique 
offers a great help to disentangle between different ex- 
tents of the polarizing sources. It is not on the scope of 
this paper to perform a detailed and complex modelling 
of several examples of sgB[c] stars. We show here the ex- 
pected impact on the SPIN signal of a typical B[e] stars 
environment, i.e. an extended, low density scattering re- 
gion, optically thin and somewhat flattened by using a toy 
model adapted from Melgarejo et al. 2001. 
For that purpose, we use the same model as for the Be 
stars with a diluted environment (typically ^ 10^ 
cm""^), more open(6' 5° — 20°) and extended (outer 
radius i?outcr ^ 40i?* — 300i?*), without any dust. The 
result for an edge-on inclination is shown in Fig. 1131 The 
principal characteristic of this environment is the angular 
diameter contrast between two distinct flux sources: the 
star and a faint but highly polarized extended envelope. 
The envelope contribution can be easily seen in Fig [TSl 
When the baseline and the polarizer are oriented both in 
the disk direction, the visibility curve is completely dom- 
inated by the stellar flux. When the polarizer is aligned 
with the p direction, the bulk of the envelope becomes visi- 



ble, superimposed on the stellar component. The envelope 
is highly resolved so that its contribution is detectable at 
a low spatial frequency: AVmax = 0.04 at /max — 0.036 
for the model presented in Fig. ^] The integrated polar- 
ization P amounts to 1.8%. The signal amplitude is close 
to one of the Be star, but its shape is very different due 
to the contrast between the point-like central source and 
the very extended and diluted environment. 
In conclusion, for B[e], we expect that the SPIN signal 
from a polarized (by electrons or dust scattering) detached 
envelope will be easier to extract from the stellar compo- 
nent due to the contrast between star and envelope ex- 
tents. The envelope polarized and natural relative con- 
tribution, its extent and geometry could therefore be re- 
trieved from a relatively simple model of its extracted vis- 
ibility curve. For Be stars such a simple reduction process 
is complicated by the fact that the disk extent is much 
lower, and because the disk is supposed to be highly op- 
tically thick for most of the models with small aperture. 

6. Instrumental Application 

After a long development, optical interferometry is now 
ready to play a significant role in astronomy. The emer- 
gence of well-funded interferometer facilities allows to 
enlarge the field of applications of this technique con- 
siderably, by increasing its reliability and its sensitivity. 
Observing in polarized light with interferometers promises 
fascinating new insights into many areas of astrophysics, 
although this capability is difficult to implement with 
current interferometers. The instrumental polarization in 
interferometers has been studied by Rousselet-Perraut 
(1996) and Elias (2001), and we suppose in this section 
that the SPIN instrument can control and calibrate the 
effect of the internal polarization on the SPIN observables. 
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Fig. 13. Model for a B[e] star environment without dust: 
A^e = 6 X 10^ (constant density), 6 = 10°, i = 90°, 

i?* = 7OR0, i?outer = 40R*, d = 2kpc. To = 0.8 in the disk 
plane, and the system is seen edge-on. We show a close-up 
of the two visibility curves \ Vs\ (solid line) and \ V-p\ (dotted 
line) . The baseline is aligned to the disk. The contribution 
of polarized light from the envelope is clearly visible in the 
oscillation of the polarized differential visibility curve AVp 
and the signal is well above the 1% sensitivity of a state- 
of-art interferometer illustrated by the dashed curve. The 
diffused light represents 8% of the total flux and is almost 
completely linearly polarized which explains why the en- 
velope signal is concentrated only in one polarization. 

In principle there is no obstacle to equip the focal in- 
strument of an interferometer by polarimetric optics. A 
simple polarization analyzer using a WoUaston prism can 
perfectly match the specifications for calibration and visi- 
bility determination in the different polarizations. In prac- 
tice we record fringe patterns in linear polarizations and 
estimate complex visibilities for each of them. The most 
promising way to calibrate the signal is to use a differ- 
ential technique (i.e. by cross-correlating the signal from 
both polarization directions for instance). The polarized 
deviation curve can therefore be determined with greater 
accuracy. It is necessary to use an unresolved star as ref- 
erence to estimate the absolute visibility. This step can be 
perfectly carried out only in natural light if the polarized 
signal can be considered as second-order effect. The first 
observations have to use large spectral bandwidth to op- 
timize the sensitivity for which a 1% accuracy is expected 
routinely in natural light. 

In the current state-of-the-art of interferometers, it is 
time consuming to record many visibility points with dif- 
ferent baseline lengthes and directions. For spherical tar- 
gets, the visibility in natural and polarized light does not 
depend on the projected baseline direction, but only on 



its length. It means that the data recorded with similar 
baseline lengths can be added in order to increase the 
SPIN SNR, since the polarizer direction of analysis follows 
the baseline movement during observation. For Be and 
B[e] stars, the measured visibility (in natural or polarized 
light) does in general not longer depend on the direction 
of the projected baseline (except for pole-on configura- 
tions) . The visibility changes with the baseline movement 
(earth rotation) which restricts the number of visibility 
points recordable per independent configuration compared 
to the spherical case. This problem is compensated by the 
large signal expected from these stars as seen in Sect El 
It can be somewhat difficult to overcome the degeneracy 
between disk density, aperture and inclination but the po- 
larized visibility provides an complementary information 
useful to constrain the parameter space. 

Which is the optimum wavelength region for this 

study? From the polarimetric point of view, the polarized 

flux generated by Thomson scattering decreases generally 

in the IR domain due to the competing influence from 

free- free continuum optical depth Tff. Since ks cx a^X^ , 

,2 

the local polarization decreases as e^"''' . For instance, 
the integrated polarization of P Cygni is decreased by a 
factor 3 between 0.55 /im and 1 /im (Nordsieck et al. 2001). 
It is therefore more interesting to observe towards shorter 
wavelengths. For Be stars the bound-free opacities cannot 
be neglected. In Tab|2 we present semi-quantitative sig- 
nal expectations for Be stars and P Cygni computed by 
means of the spectro-polarimetric data available for these 
stars. As an example, the spatially integrated polariza- 
tion of 7 Cas declines from 0.6% at the Balmer jump, to 
reach 0.52% at 0.5/im, and 0.4% at 0.66 /im. This evo- 
lution reflects only the changes in the free-bound opacity 
and free-free emission towards the envelope and affects the 
local polarization and the SPIN signal. The differences be- 
tween the spectrophotometry of 77 Tau (nearly pole-one), 
7 Cas (i~ 45°) and C Tau (nearly edge-on) are mainly due 
to inclination effects. 

The IR domain is also less attractive in term of spatial 
resolution. This is particularly striking in the context of 
hot stars since even for the examples presented in this pa- 
per, the minimum baseline range needed to resolve the po- 
larized environment is 150-200m, i.e. at the upper limit of 
possible VLTI baselines. The number of available baselines 
for such a scientific task is thus dramatically decreased. 

On the other hand, the disturbing effect from the at- 
mosphere is more striking in the optical, and the interfer- 
ometers able to perform observations in the optical wave- 
length range are currently few, and somewhat less sensi- 
tive than NIR ones. We nevertheless estimate that the gain 
in spatial resolution and polarized signal is such that the 
visible wavelength range should be preferred for SPIN ob- 
servations of hot stars environments, at least for electron 
scattering studies. We also want to mention the numerous 
applications of SPIN for the study of dusty environments 
studies, which are not directly in the scope of this paper. 
For these environments, the requirements in terms of spa- 
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Table 2. Comparison of the SPIN signal between two wavelengths. Be stars radii are estimated from Quirrenbach et 
al. 1997, and P Cygni radius from this study. The wavelength correction factors are estimated from spectropolarimetric 

observations: Quirrenbach ct al. 1997 for Be stars, and Nordsieck ct al. 2001 for P Cygni. In the /max columns the 
corresponding baselength are indicated in units of meters. In the AV50 columns the expected AVp for a 50 m baseline 
are reported. 



Name 


Distance 




Aeff 


=0.55/jm 




Aeff = 


=0.66/im 




in pc 


in mas 


^ Vinax 


/"max 




AVso 


A Vinax 


/"max 


at/so 


7 Ceis 


190 


0.56 


0.05-0.06 


100 




~0.02 


0.03-0.04 


120 


~0.015 


C Tau 


130 


0.4 


0.02-0.03 


130 




-0.015 


0.015-0.025 


160 


<0.01 


ri Tau 


110 


0.71 


0.03-0.04 


40 
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tial resolution are strongly decreased, and the IR domain 
is best suited. 

At that moment, only one long-baseline interferom- 
eter, the GET, is equipped with a polarimeter device 
for routine observations. The GI2T-REGAIN spectro- 
interferometer is composed of two 1.5-m telescopes which 
can be displaced on a North-South baseline spanning from 
12m to 65m. The REGAIN beam combiner forms the fo- 
cus in visible light and is equipped with a visible spectro- 
graph in whch a polarimetric device can be inserted. For a 
complete optical scheme and the status of the instrument, 
see Mourard et al. 2000a and Mourard et al. 2002. The 
technical description of the polarimetric observation can 
be found in Rousselet-Perraut et al. 2002. 

In the near future, further instruments will provide 
polarimetric facilities. For instance, the foreseen near- 
infrared AMBER/VLTI could observe, in a first phase, 
with a polarimetric device which processes only half of 
the incoming light. This device has been implemented in 
order to control the instrumental polarization, but this 
technical constrain could become a great opportunity to 
test the SPIN technique since it allows to extract directly 
the AVp parameter. 

A concept of a polarimetric interferometer for the 
Very Large Telescope Interferometer (VLTI), VISPER 
(VLTI Imaging Spectro-PoIarimetER) has been presented 
by Vakili et al. 2001. Since great care was taken on 
instrumental polarization effects during the design and 
construction of the VLTI, polarimetry could be straight- 
forwardly implemented in the interferometric laboratory. 
The large baselines, high-order adaptive optics (in the 
optical), fringe tracking and the foreseen dual- field 
facility PRIMA can greatly enhance the sensitivity of 
an interfero-polarimeter to the measurement of small 
polarization effects on the visibility and extend the 
number of stars for which a SPIN signal can be detected. 
Finally, the SPIN technique implies that the observa- 
tions are carried out with a certain spectral resolving 
power. This resolving power could allow a study of the 
visibility through stellar spectral lines in a similar way as 
spectropolarimetry. The differential information provided 
by a comparison between the continuum and the lines 
is richer and more sensitive than the one provided by 
classical interferometric technics. It allows to retrieve the 



differential phase information, related to the evolution 

of the stellar photocentcr through the spectrum, and 
also a dynamic information on the polarized environment 
through a spectral line (see Vakili et al. 1997, Berio et al. 
1999). With PRIMA the calibration can be even carried 
out in an absolute way. As an example, the hot star 
emission lines can often be considered as unpolarized 
and offer a good opportunity to differentially calibrate 
the polarized continuum, especially for WR stars. The 
dynamically complex environment of Be stars could also 
be studied with SPIN. Poeckert & Marlborough (1978) 
have demonstrated that the Ha line is polarized and 
contains information on the disk structure and dynamics. 
This subject should be investigated in a future work. 



7. Conclusions 

Spectro-Polarimetric INterferometry is a new and com- 
plementary way to study the polarized light from hot star 
environments. The purpose of this article is to present an 
overview of the potential of the SPIN technique, limited 
to the study of Thomson scattering within the wind of 
hot stars. The main point of the SPIN concept is that 
the star appears to have a different radius depending on 
the orientation of the polarizers relative to the baseline. 
This investigation demonstrates that state-of-the-art opti- 
cal interferometers can detect the visibility difference (the 
'SPIN signal') for a broad range of stellar wind parame- 
ters. The signal is particularly large for the denser winds 
and can be detected with small baselines (50-150m). For 
stars exhibiting dense winds, the envelope is extended and 
polarized by a large amount whereas the direct, mostly 
unpolarized stellar light is damped by the scattering, and 
only barely resolved by the interferometer. The SPIN tech- 
nique also provides a great wealth of information on Be 
star disks: disk radii in natural light for different base- 
lines, the difference of visibilities between two direction of 
polarization and baseline orientations. When the spatial 
scales of the polarized environment and the central star 
are different, the extraction of the spatial information is 
simplified, as seen for B[e] stars. 

This work has also to be seen in a larger context. The 
signal modelled in this study is certainly not the largest ex- 
pected from non-extended astrophysical targets. The envi- 
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ronments described here are geometrically and physically 
similar to astronomical objects ranging from the dusty 
environment of Yoimg Stellar Objects. AGB stars to the 
complex environment of AGNs whose apparent angular di- 
ameter is equivalent to the targets presented here. Thus, 
we strongly advocate the development of interfcromctric 
devices dedicated to SPIN measurements within the frame 
of second generation VLTI instruments. 
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